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Abstract 



We consider the orthogonal polynomials on [—1, 1] with respect to the weight 

Wc {x) = h (x) (1 — x)" (1 + x)'^ Ec (x) , a,l3 > —1, 



where h is real analytic and strictly positive on [—1,1], and Sc is a step- like function: 
Sc(a;) = 1 for X e [—1,0) and Sc(.x) = c^, c > 0, for x e [0,1]. We obtain strong 
uniform asymptotics of the monic orthogonal polynomials in C, as well as first terms of 
the asymptotic expansion of the main parameters (leading coefficients of the orthonormal 
polynomials and the recurrence coefficients) as n oo. In particular, we prove for 
Wc a conjecture of A. Magnus regarding the asymptotics of the recurrence coefficients. 
The main focus is on the local analysis at the origin. We study the asymptotics of the 
ChristofFel-Darboux kernel in a neighborhood of the jump and show that the zeros of the 
orthogonal polynomials no longer exhibit clock behavior. 

For the asymptotic analysis we use the steepest descendent method of Deift and Zhou 
applied to the non-commutative Riemann-Hilbert problems characterizing the orthogonal 
polynomials. The local analysis at x = is carried out in terms of confluent hypergeo- 
metric functions. Incidentally, we establish some properties of these functions that may 
have an independent interest. 

1 Introduction and statement of results 

1.1 Introduction 

Szcgo is the founder of the modern asymptotic theory of orthogonal polynomials on the unit 
interval for weights w that satisfy the Szego condition 



For the classical Jacobi weights the asymptotic results both on and away from the interval of 
orthogonality, as well as at its endpoints, can be derived using multiple identities that these 
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orthogonal polynomials satisfy: the differential equation, the Rodrigues formula, integral 
representation, etcetera. However, in a general situation the problem is much more difficult. 
Starting from the 80's, many new asymptotic results were found for various classes of weights, 
and the breakthrough was partially motivated by the development of the tools from potential 
theory and operator theory. 

An important new technique for obtaining asymptotics for orthogonal polynomials in all 
regions of the complex plane is based on the characterization of the orthogonal polynomials 
by means of a Riemann-Hilbert problem for 2x2 matrix valued functions due to Fokas, Its, 
and Kitaev [10] , combined with the steepest descent method of Deift and Zhou, introduced 
in [7] and further developed in [21 [6l |9] , to mention a few. 

A crucial contribution to this method is |14) . where the complete asymptotic expansion 
for the orthogonal polynomials with respect to a Jacobi weight modified by a real analytic 
and strictly positive function is obtained. However, not much is known in the case when the 
weight has a jump discontinuity on the interval. So far, the only contribution is [13 , where 
the authors considered an exponential weight on M with a jump at the origin, although from 
a different perspective of asymptotics of Hankel determinants. 

Combining ideas from [TB] and [TJj, we consider polynomials that are orthogonal on a 
finite interval [—1,1] with respect to a modified Jacobi weight with a jump, namely 

Wc{x) = il-xni + x)^h{x)E,{x), xG[-l,l], (2) 

where a, /3 > —1 and h(x) is real analytic and strictly positive on [—1, 1], and He is a step-like 
function, equal to 1 on [—1, 0) and > on [0, 1]. Observe that wi is the weight considered 
in [H] . The main asymptotic difference between the polynomials orthogonal with respect to 
wi and Wc, (c 7^ 1), lies in their behavior near the origin. While in both cases the analysis near 
the endpoints of the interval typically involves Bessel functions, only for c 7^ 1 do confluent 
hypergeometric functions appear around the origin. 

We use Pn{x) = Pn{x;wc) to denote the monic polynomial of degree n orthogonal with 
respect to the weight Wc on [—1,1], 

J Pn{x; Wc)x^Wc{x) dx = 0, for A; = 0, 1, . . . , n — 1, 

and Pn{x) = Pn{x; Wc) to denote the corresponding orthonormal polynomials, 

Pn{x^ = knPn{x), 

where A;„ > is the leading coefficient oi pn- 

The leading term of the asymptotics of polynomials p„ and Pn{z) (as n 00) for a weight 
satisfying the Szego condition ([T| (and Wc does) and z G C \ [— 1, 1] is well known, see j21| . 
It can be formulated in terms of two functions that will play a relevant role in what follows, 
and that we introduce here. Namely, 

ip{z) = z + ^Jz"^ - 1 (3) 
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is the conformal map from C \ [— 1, 1] onto the exterior of the unit circle, with the branch of 
\/ — \ that is analytic in C \ [—1, 1] and behaves like z as z — > oo. Furthermore, since the 
weight Wc on [—1, 1] satisfies ([T]), we can define the so-called Szego function -D(-z) = 
associated with tWc, given by 

D{z) = exp j ^ -^/TT^—x j ' - ^ C \ [-1, 1], 

again with ^J z"^ — \ > for z > 1 and \/l — > on (—1,1). The function D{z) is a 
non-zero analytic function on C \ [—1,1] such that 

D+{x)D-{x) = Wcix), for a.e. x £ ( — 1, 1), 

where D^(x) and D-(x) denote the limiting values of D(z) as z approaches x from above 
and below, respectively. In particular, by ([T]), the limit 

Doo = lim D{z) = exp ( — f — g^c( ) 



X- 



exists and is a positive real number. From Szego's theory (see [2T]) it follows that 



^{zY I?(z;u;c) ^/2(z2 _ 1)1/4 



[1 -I- o(l)] , as n ^ oo, (4) 



uniformly on compact subsets of C \ [—1, 1]. Using the multiplicative property of the Szego 
function, we conclude that in comparison with the case c = 1, for c 7^ 1 there is an extra 
factor, corresponding to the Szego function of the pure jump He. 

In this paper, we give uniform and more precise asymptotic results for the special weights 
([2]). We obtain the first terms of the asymptotic expansions for kn-, Pm and pn, as well as for 
the coefficients a„ and bn in the three-term recurrence relation 

Pn+l{z) = {Z- bn)Pn{z) - alPn-l{z) , (5) 

satisfied by the monic orthogonal polynomials. 

From our analysis we are also able to derive strong asymptotics for the orthogonal poly- 
nomials in the open interval (—1, 1), near the endpoints ±1, and what is most interesting, in 
a neighborhood of the origin where the jump of the weight takes place. 

Since the behavior of the polynomials P„ away from the origin is very similar to the case 
c = 1 treated in ^l4j, we will not present all formulas here. However, all the ingredients 
are contained in the results of the steepest descent analysis performed in Section [2] so that 
an interested reader can effortlessly derive the omitted asymptotic formulas. In this paper 
we concentrate on the features of the polynomials and their coefficients that stem from the 
discontinuity of the weight at the origin. 
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1.2 Asymptotics away from the interval of orthogonality 

In order to formulate our results we need to introduce some notation. For h{x) real analytic 
and strictly positive on [—1,1], n E N and c > we define the following real- valued function 
and real quantities: 

dcf Vl - x2 log hit) dt ^ I ^ ^\ (a\ 

Kx) = — ^ f , , xE(-l,l), (6) 

, , dcf loec , , , , nvr B — a ^ , , 
r]n = Vnic) = log(4n) + — + vr + h{0) , (7) 

where J is the integral understood in terms of its principal value. In general, we assume 
always Vl — 2;^ > for x G (—1,1), unless stated otherwise. 

We also introduce what will play the role of the main phase shift in all asymptotic for- 
mulas, 

f«.(.)S'2(,„-a.g(r(i^))). ifc^l, 
l2{i)„ + f), ifc=l, 

where r(-) is the Gamma function; for purely imaginary values of A 7^ 0, we take arg(r(A)) G 
(-7r/2,7r/2). 

The simplest asymptotic result concerns the monic orthogonal polynomials P„. Observe 
that a full asymptotic expansion for the usual Jacobi polynomials (/i = 1, c = 1) can be 
found in |2H Theorem 8.21.9], while for general real analytic and positive h (but with c = 1) 
it was established in |14j . Here we find only the first two terms of the asymptotic expansion, 
improving (|4]): 



Theorem 1 We have that 

2"Pn(^) ^ Poo ^{zyl^ 



as n ^ 00, 



uniformly on compact subsets o/C\ [— 1, 1]. The function 7i{z) is analytic on C\ [— 1, 1], and 
given by 

with On defined in (|8|. 

Remark 2 A more detailed analysis of the Szego function D{-] Wc) is carried out in Section 



2.3 We can simplify notation in the formula above observing that 

(Z^ - 1)1/4 
^(z)l/2 
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is the Szego function for the weig ht ^f\ — on [— 1, 1], and it takes the value 2 ^1"^ at infinity. 
Hence, 

D{z;Wc)V2{z^-iy/^ D{z-Wc)' 
where Wc{x) = \/l — Wc{x) is known as the trigonometric weight associated to Wc- 

The RH analysis performed below for z ^ [— li 1] allows also to establish a result for 
some relevant parameters associated with the orthogonal polynomials. Recall that the monic 
polynomials P„ satisfy the three term recurrence relation (|5]). The asymptotic behavior of 
these recurrence coefficients (as n — > oo) is given in the following theorem: 



Theorem ^ As n ^ oo, 



^ ^°^%in{9n) + o(^], (10) 



2 27rn \tt? 

hn = -^C0s{en)+0{\], (11) 



vrn \ n 



with On defined in 

Remark 4 In [T7j, A. Magnus studied weights of the form 

{I - x)"" [l + xf \xo - x[< X 



B, for X G [— 1, xq) 
A, for X G [xo, 1] , 



with A and i? > and a, j3 and 7 > —1, and xq G ( — 1,1). Formulas (lO)-(ll) show that 
for 7 = Xq = the asymptotic behavior of the recurrence coefficients conjectured in |T7| is 
correct, with the possibility to replace o(l/n) by 0(l/n^) in the error term. For more details 
see Section 3.2 below; the proof of the conjecture in its full generality is contained in fll]. 

The leading coefficients kn of the orthonormal polynomials pn satisfy the following asymp- 
totic relation: 



Theorem h As n ^ 00, 



nn 

k 



with On defined in (|8| 



2a2 + 2/3^-1 log(c) /log(c) ^^^ 1 n n 



5 



1.3 Local asymptotics 

Now we need to introduce further notation. Set 



G(a;C) = iFi(a;l;C) e 



-C/2 



(A;!)2 ^ 



E 

k=0 



(12) 



where iFi (a; 6; •) is the confluent hypergeometric function; G is an entire function of C for 
any value of the parameter a S C, and G{a; 0) = 1. Furthermore, for x E {—5, 0) U (0, (5) let 



p{x) 



log c 



TT 



log 



arcsmlx 



2x 



1 — x"^ 



a + [3 



arcsin(a;) + h[x) — h{0), 



(13) 



(14) 



completed to a continuous function on (—1, 1) by p(0) = 0. 
Set also 

T(c)^=^sgn(log(c))y^^, c^l, T(l) = l, 

where we always take the positive value of the square root. The asymptotic behavior of P„ 
on compact subsets of an interval {—5, 6) C (—1,1) is given by the following theorem: 

Theorem 6 For 6 £ (0,1), locally uniformly on compact subsets of {—5,6) the following 
asymptotic formula holds: 



Pn{x) 



T(c) 



171—1 



/2y^^i;^(l-x2)l/4 



X Re 



-7r~arcsin 



(A; 2in arcsin (x)) ( 1 + 



Unix) 



+ 



n 



with 



Tin {x) 



4a2 - 1 



8(e^ 



+ 



4/32 



(3;) _ 1) 8(^gjarccos(a:) _|_ ]^-) 



log(c) 



27rX6* 3.rccos(, 



( cos(6'„)e 



i arccos(a;) 



+ sin(6'„) 



log(c) 



TT 



+ 



i log c / log c 



+ e 



27rarcsin(x) \ vr 
p{x) given in (13), 6'„ in ([s]), and A = ilog(c)/7r. 

Remark 7 The Riemann-Hilbert analysis we perform next gives us an asymptotic expression 
for PnS in a small disk of the complex plane centered at the origin, see formula ( 101 ) in Section 
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Corollary 8 Locally uniformly for x G (—5, 6), 6 G (0, 1), 



\n J 2"-i/2ycr^o) 
^„ given in ([8| and A = ilog(c)/7r 



Im 



(15) 



See Figure [T] for a typical behavior of the function in the right hand side of (15) close to 
the origin. 




Figure 1: Typical graphics of the r.h.s. of (15) near the origin. 



Recall that Pn has n simple zeros, all lying on (—1, 1). It is well known that they distribute 
asymptotically in the weak-* sense according to the equilibrium measure of the interval. In 
other words, the normalized zero counting measure for the sequence Pn weakly tends to the 
absolutely continuous measure on [—1, 1] given by uj(x) dx, with 



a;(x) 



of 1 



x^ 



As it follows from several works of Deift and collaborators (and also from a recent series of 
papers of Lubinsky and Levin and Lubinsky, see e.g. |15| I16j). a much stronger statement 
holds: at any point of (—1,0) U (0, 1) they distribute very precisely in accordance with u!{x), 
complying with the so-called "clock behavior", see e.g. [19]. If, following [[19j, we enumerate 
the zeros x^"^ of Pn as follows, 

•••<x(1<...< 



< xLI < < xj,") < 



< X 



in) 



< . 



(16) 
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then "clock behavior" at the origin (where i^(0) = l/vr) means 

hm - (xgi-x;.")) =1, jeZ. (17) 

Proposition 9 If c > 1, then the sequence {nxQ^^/vr} is dense in an interval of the form 
[0,t], where t = t{c) < 1. Furthermore, 

< hminf - (xP - x^f^}A < hmsup - (x^!"^ - < 1, ken, 



and 



> 1, -k£N. 



n TT 

In particular, the clock behavior of the zeros of Pn at the origin does not hold. 
IfO<c<l, the same inequalities hold inverting the roles of k and —k. 

This result is not surprising, taking into account that x = is not even a Lebesgue point for 
the weight Wc, that is, regardless of the meaning we give to Wc(0), 



lim 



1 

- / \wc{x) — Wc{^)\ dx ^ ^. 

S J ~s 



However, to the best of our knowledge, voc with c 7^ 1 provides the first instance of an explicit 
orthogonality measure for which the clock behavior fails in the bulk (interior of its support). 



Remark 10 A weaker condition than (17 1 is the quasi-clock behavior (see [IH])) namely 

(n) (n) 

This limit is violated in our situation too. However, 

lim liminf— (x^^2^ — x^f^] = lim limsup— \ x^^2^ — x^f^] =1, (18) 

which shows a smooth transition to the genuine clock behavior as we move away from the 
jump of the weight. 

Very much related with the clock behavior is the "universality problem" for the Christoffel- 
Darboux (or CD) kernel 

n-l 

Kn{x,y) = '^Pk{x)pk{y), (19) 

fc=0 

where pn are the orthonormal polynomials with respect to the weight Wc- This problem 
has its origin in the random matrix theory and has been attracting lately close attention of 
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many researchers. A recent series of remarkable contributions of Lubinsky allowed to weaken 
considerably the conditions on the weight to be able to assure universality: now we know 
that for t within the support of the weight where it is continuous, 



lim — 1^=^= 



TTX 



,t + 



Try 



n 



n 



sin {^{x — y)) 
^{x - y) 



(20) 



The right hand side is the well-known sine (or sine) kernel; for our weight Wc, this formula is 
vahd for t G (-1, 0) U (0, 1). It was observed in [E] that (|20j) imphes ([T7|. 

We show that the jump discontinuity in the weight leads to a different kernel, constructed 
in terms of the confluent hypergeometric function defined in (12): 

Theorem 11 For c > 0, c 7^ 1, locally uniformly for x and y on {—5,5), < 5 < 1, 



IT fiTX -Kyx 
hm - A„ — , — = Aoo [x, y) , 
\ n n J 



(211 



with 
Koc {x, y) 



1 logc [G (1 + A; 27rix) ; G (A; 27riy) " 



h{0) nic^ -1 
2 log c 



X -y 



X j^y, 



(22) 



h{0) C2 



1 



(G' (1 + A; 27rix) G (A; lirix) - G (1 + A; lirix) G' (A; 27rix)) , x = y, 



where A = ilog(c)/7r, G was introduced in (12), and as usual, [f{x);g{y)] = f{x)g{y) 

f{y)g{x)- 



Several remarks are in order. 



Since G' (1 + A; 0) = A + 1/2 and G (A; 0) = 1, evaluating Koo{0, 0) in (22 ) we conclude 



that 



lim 

n-.^ 0,0 



TX \ 

n J 



G' (1 + A; 2TTix) G (A; lirix) - G (1 + A; lirix) G' (A; Inix) , 



locally uniformly in {—5,5). This shows that even the weak Lubinsky's "wiggle condition" 
(term coined by B. Simon, see e.g. |19| Theorem 3.6]) is not satisfied in a neighborhood of 
the jump of the weight. 

The kernel for x 7^ y in (22) is written in the so-called integrable form. Taking into 
account the properties of the functions in the right hand side, we can rewrite it alternatively 
in a totally real form: 



Koo{x,y) 



log c 



Tr{x - y) h{0) - 1 



lm{G{l + X;2TTix)G{X;2TTiy)) , x + y. (23) 
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Since G{l,z) = exp(z/2), straightforward computations show that as c — > 1, K^o reduces 



to the sine kernel. Notice that combining ideas from |16! and [18] we can use (23) to arrive 
at the same conclusions about the spacing of zeros of P^s as we did at the end of Subsection 

El 

The confluent hypergeometric functions appeared in the scaling limit (as the number of 
particles goes to infinity) of the correlation functions of the pseudo-Jacobi ensemble in [3]. 
This ensemble corresponds to a sequence of weights of the form 

(1 + x2)-"-^<^We2i""(")'''^s(^+^^\ x£R, (24) 

where n is the degree of the polynomial and s is a complex parameter. The connection between 



both problems becomes apparent if we perform the inversion x ^ 1/x in ( 24 ) ; this creates at 
the origin an algebraic singularity with the exponent Re(s) and a jump depending on Im(s). 
Koo is a particular case of the reproducing kernel obtained by Borodin and Olshansky in 
Theorem 2.1 of [Sj when Re(s) = 0; for a general situation, see ^llj . 

A recent paper of Lubinsky [15] revealed an interesting connection of K^o with the theory 
of entire functions. Namely, in accordance with Theorem 1.6 of [15j, K^o is a reproducing 
kernel of a de Brange space, equivalent to a classical Paley- Wiener space. More precisely and 
following the notation of [15j, the Hermite-Biehler class HB is the set of entire functions 
E with no zeros in the upper half plane C"*" = {Imz > 0} and such that |-E(2;)| > |-£'(^)| 
for z G C^. The de Branges space 'H{E) corresponding to -B G HB is comprised of entire 
functions g such that both g(z)/E{z) and g{z)/E(z) belong to the Hardy class H'^{C'^). A 
reproducing kernel for TCi^E) is 



^ i E{x)E{y) - E{x)E{y) 
/C(x, y) = ^ , x^y. (25) 



Comparing this expression with K^o in (|22| we conclude that 

/C(x,y) = i^oo(a;,y). 



with A = ilog(c)/7r and 



log c 



1/2 



(see below). Lubinsky showed that reproducing kernels, different from the right hand side in 



(20), can appear for sequences of measures (cf. p]). To the best of our knowledge, this is the 
first explicit example of a non-sine reproducing kernel of a de Brange space that arises as a 
universality limit in the bulk of a fixed measure of orthogonality. 

The proof of the asymptotic results stated in this paper (see Section [3| is based on the 
steepest descent analysis of the Riemann-Hilbert problem that we carry out in Section [2j 
A key step is the construction of the local representation at the origin, which is done in 
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Subsection 2.6 The study of the zeros of P„'s at the origin, the analysis of the clock behavior 
and the connection with the de Brange spaces requires some further properties of the confluent 
hypergeometric function (A; 1; 2;), which we were unable to find in the literature and which 
might have an independent interest. We summarize them in the next proposition; the proofs 
are relegated to Section |4j 

Proposition 12 Let a£R\ {0}. Then 
(i) functions 



fi{z) = G{ia,iz) and f2{z) = G{1 + ia,iz) 



(see (12)) belong to the Hermite-Biehler class HB; 
(ii) for X gM, 

iFi(m;l;ix) /O and iFi (1 + ia; 1; ix) / 0. (26) 

In particular, all zeros of iFi {ia; l;iz) lie in the lower half plane C~ = {Imz < 0}, 
while the zeros of iFi (1 + ia; 1; iz) lie in the upper half plane C"*". Additionally, 

\iFi{l +ia;l;iz)\ < \iFi{ia;l;iz)\, Imz > 0, (27) 

and the equality holds only for z gM. 

(Hi) if a > 0, the function 

y(x) = argiFi (ia; l;ia;) , y(0) = 0, 

is real- analytic and non-positive, strictly increasing on the negative and strictly decreas- 
ing on the positive semiaxis. For a < the same assertion is valid replacing y{x) by 
—y{x) . It is also the solution of the following initial value problem: 

xy' = a (cos (x - 2y) - 1) , 2/(0) = 0; (28) 

(iv) for a G M, the function 

(J5(x) =?x-2 arg(iFi(m;l;ix)) = x-2y(x), 6(0) = 0, (29) 
is strictly increasing in M. 



Remark 13 The assertion in (i) does not imply that iFi {ia; l;iz) G HB, and in general, 
this is not true. 

Interestingly enough, the proof of (i) is based on some properties of the Christoffel- 
Darboux kernel observed by Lubinsky in [15]. In this sense, the theory of the confluent 
hypergeometric functions has benefited from the properties of the reproducing kernels. In 
the opposite direction, the strict inequality in (27) implies that Koo{z,z) > for z G C \ M 
and c / 1, see (113) below. 
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2 The steepest descent analysis 



2.1 The Riemann-Hilbert problem 

Following Fokas, Its and Kitaev ^lOj we characterize both the orthogonal polynomials and the 
CD kernel in terms of the unique solution Y of the following 2x2 matrix valued Riemann- 
Hilbert (RH) problem: for n G N, 

(Yl) Y is analytic in C \ [-1, 1]. 

(Y2) On (—1,0) U (0,1), Y possesses continuous boundary values Y+ (from the upper half 
plane) and Y_ (from the lower half plane) , and 

Y,(x)=Y_(.)(; -f)). 



(Y3) As z ^ oo. 



where I is the identity 2x2 matrix. 



Y{z) ={l + 





z^" ' ' 



(Y4) Y has the following asymptotic behavior at the end points of the interval: for G 
{—1, 1} set s = a if C = 1, and s = /3 if C = —1- Then for z ^ z £ C\ [—1, 1], 



Y{z) = { 



O 



1 \z-cn 

1 k-CIV' 



if s < 0; 



Ol' , if. = 0; 

1 iog|z-cir 



1 1 
1 1 



if s > 0. 



Furthermore, at the origin Y has the following behavior: for z ^ 0, z € C\ [—1, 1], 

\1 log\z\ 



Standard arguments (see e.g. [E]) show that this RH problem has a unique solution given by 



Y{z,n) 



PrAz) 



C (PnWc) {Z) 



-27rifc2_lP„_l {z) -2^ikl_^C {Pn^lWc) {z) 



(30) 
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where P„ is the monic orthogonal polynomial of degree n with respect to Wc', kn is the leading 
coefficient of the orthonormal polynomial pn, and C (•) is the Cauchy transform on [—1,1] 
defined by 

Clearly, Y and other matrices introduced hereafter depend on n, fact that we indicate writing 
Y(-,n). However, we omit the explicit reference to n from the notation whenever it cannot 
lead us into confusion. 



2.2 First transformations 

We apply the Deift-Zhou method of steepest descent to the RH problem above; some of the 
steps are standard and we occasionally omit those less relevant details that can be easily 
found in literature (each time we try to provide a suitable reference though). As in ([s]), 
ip denotes the conformal mapping from C \ [—1,1] onto the exterior of the unit circle. Let 

(73 = be the third Pauli matrix; in what follows, for a G C\{0} we use the notation 



dc! I a 
1/a 



then for b G C, a^^3 is understood as (a^)*^^. Furthermore, if 7 is an oriented Jordan arc, and 
an analytic function / has boundary values at 7, we denote by /+ (resp., /_) its boundary 
values on 7 from the left (resp., from the right). 
Set 

T (z) = 2'"'^Y (z) if {z)-'"'\ (31) 
Then T is the unique solution of the following equivalent RH problem: 

(Tl) T is analytic in C \ [-1, 1]. 

(T2) On (—1, 0) U (0, 1), oriented from — 1 to 1, T possesses continuous boundary values T_|_ 
and T_, and 

V+^" (x) Wc (x) 

99?" (x) 



T+(x) = T_(x)j ^ 



(T3) As z ^ 00, 

T{z) =1 + 



(T4) T has the same asymptotic behavior as Y at ±1 and . 
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Next transformation is based upon the factorization of the jump matrix for T: 



In order to introduce a contour deformation we need to extend the definition of the weight 
of orthogonahty to a neighborhood of the interval [—1,1]. 

By assumptions, /i is a holomorphic function in a neighborhood U of [—1, 1], and positive 
on this interval. We set 

w{z) = h{z) (1 - (1 + zf (33) 

holomorphic in U \ ((— cxd,— 1] U [l,+oo)), and such that w{x) > for x E ( — 1,1). In 
particular, w(0) = h(0). We also extend the definition of the step function He by 




if Re z < 
if Re z > 0. 



Then we set 

w,{z)'^w{z)E,{z), (34) 

which is a holomorphic function in U = U \ ((— cxd, —1] U [1, +oo) U iR). 

With this definition the left and rightmost matrices in ( |32[ ) have an analytic extension to 
the portion of U in the lower and upper half plane, respectively, and we can define the next 
step: lens opening or contour deformation. Namely, we build four new contours 7i lying in 




Figure 2: First lens opening. 

U (except for their end points) such that 71 and 73 are in the upper half plane, and 71 and 
72 are in the left half plane, and oriented "from —1 to 1" (see Fig. [2|. This construction 
defines three domains: the inner upper domain, bounded by [—1, 1] and the curves 71 and 73; 
the inner lower domain, bounded by [—1, 1] and the curves 72 and 74, and finally the outer 
domain, bounded by curves 7^ and containing the infinity. 
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Using the matrix T from (31 ) we define the new matrix S by 

'T (z) , for z in the outer domain, 

0^ 



S (z) = <^ 



T(z) 



1 0\ 



for z in the inner upper domain, 
for z in the inner lower domain. 



Then S is the unique solution of the following RH problem: 

(51) S is analytic in C\S, where S = [—1, 1] U IJi=i 7*- 

(52) S satisfies the following jump relations: 



S+(a;) = S_(x) 



1 

1 / \-2n 



Wc ix) 
1 



, for j^Uj^j 

o""^), for xG (-1,0) U (0,1). 



(S3) As z ^ oo. 



S(z) =i + o 



(S4) S has the following asymptotic behavior at the end points of the interval: for C ^ {" 
set s = a if C = 1, and s = (5 \i Q = —\. Then for z ^ C, z G C \ S, 



for s < 0: 



for s = 0: 



for s > 0: 



S(z) = O 
S(z) = O 



1 i^-cr 
1 k-cr 



as z — > C; 



log|z-C| logjz-Cl 
log|z-C| logjz-Cl 



as z ^ C; 



S(z) = <^ 



o 
o 



as z ^ C from the outer domain; 
, as z ^ C from the inner domains. 
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(S5) S has the fohowing behavior at the origin: as z ^ 0, z G C\S, 

( 

I I Incr \7.\ \ 

as z ^ from the outer domain; 



S(z) 



O 

o 



1 log \z\ \ 
1 log \z\ I ' 
'logl^l log|z| 
log|z| log|z| 



as z ^ from the inner domains. 



2.3 The Szego function for Wc 

In this section we analyze in detail the structure and properties of the Szego function in- 
troduced in Section Recall that for a non-negative function h on (—1,1) satisfying the 
Szego condition 

^1 logh{t) , 

we define in C \ [—1, 1] its Szego function D(-, h) by 



D{z, h) = exp 



V^23T /■! \ogh{t) dt 

2vr ~t 



exp 



loghjt) 



(36) 



with (VT" 



> for z S ( — 1, 1) in the rightmost expression in (36). 



Due to the multiplicative property of the Szego function, we have that for Wc defined in 

D{z,Wc) = D{z,w)D{z,Ec) . 



(37) 



Straightforward computation shows that 

(^_l)"/2(^ + l)/3/2 



D{z,w) = D{z,h) 



if> 2 (z) 



Z)(2;, Sc) = c exp —A log 



1 - - 1 



where D{-,h) is computed by formula (36), and 



dcf .logc 

A = I 

vr 



(38) 
(39) 



We must clarify that in (38) we take the main branches of {z — 1)"''^, {z + 1)^^"^ and — 1 
that are positive for z > 1, as well as the main branch of the logarithm. 



From (37) we obtain that 



Doo = D{oo, Wc) = \fcD (oo, K) 2 



^(q+/3)/2 



> 0. 



(40) 



Let us study the boundary behavior of the Szego function on the interval. By ( 38 ) , 

a+|3 



lim 7?(z,zi;) =e™/^(l-x)"/2(l + x)''/V, Ma;) Hm B{z,h), 

z^xe(-l,l), z^iG(-l,l), 

Im2>0 Im2>0 
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where 

^^{x) = x + iVl-x2 = e'^''""'<''\ (41) 

with (1 — x)"'^^, (1 + x)^^^ and \/l — positive for x G (— 1, 1). 
Analogously, 

lim D{z,w) =e-^''''^{l-x)'^/^ {l + xf^ ^^{x) lim i:>(z,/i). 

2^a;g{-l,l), z^a;6(-l,l), 

Im2<0 Im2<0 

We can be more specific about the limit values of D(z,h) on (—1,1) if we use the 
Sokhotskii-Plemelj formulas [,12^ Section 4.2]: 

/ log/i(t) \ l log/i(t) 1 f} \ogh{t) dt 

where J is the integral understood in terms of its principal value. So, if we define ^(x) as in 
^ , then using ( 36 ) we get 



^-ih(x) 



lim D (z, h) = \/h{x) 

lmz>0 

lim L»(z,/i) = y^e^^(^). 

z->a;e(-l,l), 

Im2;<0 

Observe that ^(x) is real-valued on (—1, 1), so that |e^''^^^)| = 1. So, if we define on (—1, 1) 
the real-valued function 

$(x) == ^-^ arccos(x) — h{x), (42) 

then 



lim D{z, w) = \/ w{x) exp {±.i^{x)) . 
z->xe(-i,i), 

±Im2>0 

On the other hand, it is easy to check that with the specified selection of the branch of 
the square root, 

1 - Wz'^ - 1 
Z 

z 

is a conformal mapping of C \ [— 1, 1] onto the lower half plane, such that the lower shore of 
(—1,1) is mapped onto itself, while the upper boundary is mapped onto (— oo,— 1) U (l,oo). 
In particular, 

lim arg = (J, lim arg = — tt. 

z^xe{0,l), \ Z I z-.a;G(-l,0), \ Z j 

ImzT^O ^ / Imz^O ^ / 



17 



Hence, 



lim D{z,r.c) = c exp —Alog 

z-.a:S(0,l). \ 

±lmz>0 ^ 



1 ± y/Y 



X 



c exp ^pAlog 



1 + vr 



with Vl — 2;^ > on (—1,1). Taking into account that e ^'^^ = c, we also get 

= exp TAlog 



lim D{z,r.c) = exp —Alog 

z^a:e{-l,0), \ 

±Im2>0 \ 

Both identities can be summarized by 



1 ± Vl -x2 



X 



1 + Vl - x^ 



X 



lim D{z, Ec) = \/eJx) exp =F 

±Im2>0 ^ 



logc 

^ log 

vr 



1 + 



In order to clarify the local behavior of D{z, He) at the origin we observe that for z G 
C\ (—00, 1] function D{z, He) coincides with 



c exp ^— Alog(l — i\/ z^ — 1) + Alog(/ 
if we take there the main branch of log(z), so that 



Since 



it yields 



-Aiog(2)^ (z, He) = c exp (-Alog (1 - _ 1 
lim log (1 — i\/ z'^ — 1^ = log(2) , 



Im2>0 



D (z. He) = ci+- (1 + (1)) , as z ^ 0, Im z > 0. 



The case Im z < can be deduced using the symmetry of D (•, We) with respect to 
We can summarize our findings in the following lemma: 



Lemma 14 The Szego function D{-, w) for the weight w defined in ( 33 ) exhibits the following 
boundary behavior: 

lim D{z,w) = \J w{x) exp (ibi<I>(x)) , (43) 



2^a:e(-l,l), 

±Im2>0 



with the notation introduced in ^ and (42). 
Furthermore, for the step function He, 



lim 

z^a:e(-1.0)U(0,l), 

±Im2>0 



i:)(z,He) = \/He(a;) 



exp =Fz 



log c 



vr 



log 



1 + ^/\~^ 



and 



D {z, He) = c^^i (1 ^ ^ (1)^ ^ asz^O, ±lmz>0. 



(44) 
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Obviously, the boundary behavior of the Szego function D{-,Wc) at (—1, 1) can be deduced 
from this Lemma and (37). 



2.4 Outer parametrix 



Since \(f {z)\ > 1 for z G C\ [—1, 1], the matrix S introduced at the end of Subsection 2.2 has 
jumps across each contour 7^ that are exponentially close to I, as long as we stay away from 
the singularities ±1 and 0. So, we can expect S to behave similarly to the 2x2 solution N 
of the following RH problem in this region. 

(Nl) N is analytic in C\ [-1, 1]; 

(N2) N satisfies the following jump relations on (—1, 0) U (0, 1): 

Wc (x) 



(N3) As ^ ^ 00, 

N(z) =1 + 



^z ^ 

An explicit solution of this problem is well-known (see e.g. jS] and |14| Section 5]) and can 
be built in terms of the Szego function D(-, Wc) and its value at infinity defined in (40): 

N{z)'^D2A{z)D{z,w,r^\ (45) 

where 

and we take the principal branches in such a way that a is analytic in C\ [—1,1] with a (2) — > 1 
as z — 00. For future reference it is convenient to notice that an alternative expression for 
the entries of A can be obtained using that 

a(z) +a-i(z) 93(2)^^^ 
A\\[Z — — 



2 ^/2fz2_ 1)1/4' 



a (z) — a 1 (z) ii^ (z 



,-1/2 



^ ' 2i V2(z2-l)V4 ^{z) 

where we take again the main branches of the roots. 

It is known (see [Ijj) that N does not match the behavior of S at the endpoints of the 
interval [—1,1], requiring a separate analysis there. Moreover, comparing the local condition 
(S5) for S with the behavior of D(-,Wc) at the origin (see (44)) we conclude that a local 
analysis will be needed also at z = 0. 
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2.5 Local parametrices at the endpoints of the interval 

We fix a (5 € (0, 1/8) and for each C G { — 1, 1} we consider the neighborhood = {z £ C : \z — (\ < 6} 
such that lies entirely in the domain U of analyticity of h. We construct a 2 x 2 matrix- 
valued function in [/^ \ S that exhibits the same jumps on T, nU(^ and the same local 
behavior at z = C as S, and that matches the matrix N on the boundary dU(^. Namely, 

(P^l) P^; is holomorphic in C/(\S and continuous up to the boundary. 

(P(^2) P^ satisfies the following jump relations: 



(P^4) has the following behavior as z — > C; ^ ^ ^c\^' '^ith s = a if = 1 and s = /? if 




(P(^3) As n ^ oo, 




uniformly for z G dU(^\T,. 



c = -i 



• for s < 




• for s = 




• for s > 




as z 



as z 



( from the outer domain; 



C from the inner domain. 



We skip the details of construction of P±i, that can be found in |14j . 
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2.6 Local parametrix at the origin 

We fix a 5 G (0, 1/8) and consider the neighborhood Uq = {z £ C : \z\ < 6} such that Uq lies 
entirely in the domain U of analyticity of h. We construct a 2 x 2 matrix- valued function Pq 
in C/q \ S that exhibits the same jumps on S n C/q and the same local behavior at z = as S, 
and that matches the matrix N on the boundary dUo. Namely, 

(PqI) Pq is holomorphic in Uq\T, and continuous up to the boundary. 

(Po2) Pq satisfies the following jump relations: 

Po+(z) = Po-(^) J(^)-2n fj , for z G [/o n (\J^^i^ \ {0}; 

(0 w (x)\ 
i_ nM, for XG [/on ((-1,0) U (0,1)). 



(Po3) As n ^ cxD, 



Po(z)N"^ (z) = I + O ( ^ ) uniformly for z G dUo\^. 



(Po4) Pq has the following behavior as z — > 0, z G Uo\T,: 



1 log |z| \ 

O [ , as z — > from the outer domain; 

1 log|z|y 

O \ ^ \ , asz— >0 from the inner domain. 
log|z| log|z|/ 



We build the solution of this problem in two steps. First we obtain a matrix P^^^ that 
satisfies conditions (PqI, Po2, Po4), and after that, using an additional freedom in the con- 
struction, we take care of the matching condition (Po3). 

Let us define at this point an auxiliary function W holomorphic in U\{{—oo, —1] U [1, +cxd)) 
given by (see ( 33 ) ) 

W{z) = y^cw{z) , such that VF(x) > for X G (-1,1). (48) 

Then 

l^(x) = |v ^^^^' -l<x<0; ^^g^ 

[y/Wc (x) C^l, < X < 1. 

We construct the matrix function Pq in the following form: 

Po (z) = E„ (z) P(i) (z) W (z)-'^^ if (z)-"'^^ , (50) 
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where E„ is an analytic matrix- valued function in Uq (to be determined). Matrix P^-*^) is 
analytic in Uo\'^', using the properties of W and ip it is easy to show that 



1/ci 
-c y 



-1/c 0, 



, xe(-<5,0), 



(51) 



and 



p« {z) = P'y (z) { 



(1) 



1 0^ 



2 e (7iU72)n[/o\{o}, 

z G (73U74)nC/o\{0}. 



(52) 



Taking into account that W (z) = O (1) and ip{z) = O (1) as z — > 0, we conclude also from 
(Po4) that P^^^ has the following behavior at the origin: as z ^ 0, z G C\S, 



p(i)( 



O 



O 



1 log|2;|\ 
^1 log|z|y ' 

logjzl log|z| 
logl^l log|z| 



from the outer domain, 



from the inner domain. 



(53) 



In order to construct p(^) we solve first an auxiliary RH problem on a set Sij, == IJj=i 
of unbounded oriented straight lines converging at the origin, like in Fig. |3j More precisely. 



Ti = |te^^/2 : t > o} , Ta 
r4 = |te^*^/2 . i > o| , Ts 



These lines split the plane into 6 sectors, enumerated anti-clockwise from ® to © as in Fig.jsj 
We look for a 2 x 2 matrix valued function ^ (z), satisfying the following conditions: 

(^1) * is analytic in C\i;^. 

('J'2) for /c = 1, . . . , 6, satisfies the jump relation ^'^(C) = ^^{Q)Jk on T^, with 



c 
-1/c 



J2 — Jf, 



1 

1/c ly' 



1 

c 1 



Ja 



1/c 
-c 



(54) 
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Figure 3: Auxiliary contours S 



(^'3) ^ has the following behavior as C ^ 0: 

'log Id logicr 

"I 

*(C) 



,iogici log Id, 
'i iogid\ 
1 log Id/' 



for C G ® U 
as C £ 



U ® U I 



If we use the notation X = i log(c) /vr introduced in ( 39 ) , then we readily see the connection 
of the RH problem above with that studied recently in |T3] . Following the approach of [13] 
(with slight modifications), we construct ^ explicitly in terms of the confluent hypergeometric 
functions 



0(a,l;C) = iFi(a;l;C) and ^ (a, 1; C) = 2^0 (a, a; 



-1/C) 



that form a basis of solutions of the confluent hypergeometric equation (11!"+ {1 — C) w' —aw 
0, see [U formula (13.1.1)]. Namely, let 

G (a; C) = </> (a, 1; C) e'^'\ H (a; C) = (a, 1; C) e~^'\ 
They are solutions of the confluent equation (see e.g. [1| formula (13.1.35)]) 

1 c 



^ II \ I \ 

Qw + + 



a w 



0; 



(55) 
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multivalued, and we take its principal branch in — | < arg {() < For these values of we 



in fact, G{a; •) is the only entire solution of (55) such that G(a; 0) = 1. Function H{a,C) is 
multiv 
define 

*(C) 



r(l-A)G(A;C) -HiX;C) 
r(l + A)G(l + A;C) 9gf//(l + A;C), 



By (^'2), if we set 



*(C) JeJi, 


for C G ® 




for C G @ 


^ {0 J(,JlJ2J^\ 


for C G @ 




for C G ® 


*(0 


for C G © 




for C G ® 



*(C) 



then ^ has the jumps across specified in (^'2). Explicitly 



*(C) 
*(C) 
*(C) 
*(C) 
*(C) 
*(C) 



r(~A) 

r(l-A)G(A;C) 
r(l + A)G(l + A;C) 

cH (A; e-2^*C) 
-c^ F (1 + A; 



^fef)//(l-A;e-C)^ 
i/(-A;e--C) 



Cg®, 



r(A) 



A;e— C) ] 
-A,l;e— C) J 

^H{l-X;e~-X) 
H{-X;e-^X) 



CG 



CG 



-^(A;C) 



cH{-X;e-X) ¥^i/(l + A;C), 



CG 



r(l-A)G(A;C) 
r(l + A)G(l + A;C) 



-i r(i-A) 
r{A) 



r(-A) 
H{1- A;e-^X) 
;-A;e"-C) 



-(-A) 

-^(A;C) ^ 
^ii±^F(l + A;C) 



CG 



-i?(A;C) 



CG 



(56) 
(57) 
(58) 
(59) 
(60) 
(61) 



Direct verification shows that ^ coincides, after an appropriate change of parameters and 
a multiplication from the left by the constant matrix 

1/c 
1 

with the solution of the corresponding RH problem found in [13] (cf. formulas (7.26) and 
(7.27) therein). In consequence, the matrix-valued function ^ defined in (56)-(61) solves the 
RH problem (^'l)-(^'3); moreover, det = 1. 
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In order to construct the analytic function E„ in (50) we need to study also the asymptotic 
behavior of ^ at infinity. Let us introduce the notation 

r(-A) 



Then for purely imaginary values of A 7^ 0, 



r(A) 



m 

T{xy 



FA 



1, 



T-A = TA- 



This value is not defined for A = 0; by continuity, we set tq = —1. 

Lemma 15 As ^ 00, C £ C \ and with the notation A = ilog(c)/7r, we have 



*(C) 



1 + 



A 



-A 



l^~A(73g-Co-3/2 




ICI 

i/f <argC< f ; 
if - ^ < argC < 



(62) 



where we use the main branch of C, ^ = e •^^"S'' with the cut along iM_. 

This result is a direct consequence of formulas (4.60)-(4.63) from [.13,, and can be obtained by 
straightforward computation using the asymptotic properties of the confluent hypergeometric 
functions (see e.g. [H formulas (13.5.1-2)]). In fact, formulas in [Ij give us the complete 
expansion of ^ . 

Now we are ready to build P^^^ as in (50). Recall that is a conformal mapping from 
C \ [—1,1] onto the exterior of the unit disk, so we can define in C \ M the analytic function 

cf J TTZ — 2 log if [z) , for Im z > 0, 



f{z) 



iri + 2 log (p (z) , for Im z < 0, 



(63) 



where we take the main branch of the logarithm. Using that (x) ip- {x) = 1 on (—1,1) we 
conclude that /+ (x) = /_ (x) there, so that / is holomorphic in C \ ((— cxd, —1] U [1, +cxd)). 
For \z\ < 1 we have 

^ ■ ' " ' 0. (64) 



/ {z) = 2iz + ^ iz^ + O {z^) , as z 



Hence, for 5 > sufficiently small, / is a conformal mapping of Uq. Moreover, by (41), 

/(x) = 2i arcsin(x), xG( — 1,1), 



(65) 



so that / maps the real interval (—1,1) one-to-one onto the purely imaginary interval (— vrz, iri). 
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We can always deform our contours 7^ close to z = in such a way that 

/(7in?7o) cTg, /(72nf/o) cTs, /(73n[/o) cr2, 7(74 n t/o) c Tg. 

With this convention, set 

P«(z)^='*(n/(z)), zeUo. 



(66) 



By (\I'l)-(^'3) and (64), this matrix-valued function has the jumps and the local behavior at 
z = specified in (51|-(53). Taking into account the definition (63) we get that 



[z), for lb Im z > 0. 



Hence, by Lemma 15 



*(n/(z)) 



1 + 



nf{z) \-1/tx a 



+ 



1 



(n/(z))-^'^3i-"'^3 



X < 



\l 



if Imz > 0; 
V3"'^3(z), if Imz < 0, 



(67) 



where the main branch of [nf {z)] is taken with the cut along (— oo, 0]. Since 

[/(z)]^ = |/(z)|^expf-i^arg(/(.)) 



vr 



straightforward computations show that 



'|/(^)l^c-V2^ forO<x<l, 
|/(x)|^c-i/2^\ for - 1< x < 0, 



where we assume the natural orientation of the interval. 



(68) 



Now we will build the analytic matrix E„ in (50). In order to comply with condition 
(Po3) above, we need 



E„(z) ~N(z)v9(z)'^"«Ty(z)"^ PW( 
uniformly for z G dUs\T,. Taking into account ( |67[ ), we define 



E„ (z) = N {z) W {z^' X <^ 



r"-3 (n/(z)) 



1 

-1 0, 



if Imz > 0; 



, if Im z < 0. 



(69) 
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By construction, E„ is analytic in f7o\M. Furthermore, by (N2) and (49), for x G (—(5,0) U 
(0,5), 



W {x)-"' Nz\x)'N+{x)W {x 



1 0-3 



Wc{x)/W^{x] 
-W'^{x)/wc{x) 

c±i- 
-c^i 



for ± Rex > 0. 



From (68) and (69) it follows that 

(x) En+ (x) = I, for X G (-(5,0) U (0,6). 
So, the origin is the only possible isolated singularity of E„, in Uq- 
Proposition 16 



limE„(z) = ^DS 



1 1 
-1 1 



with rjn introduced m ([7]). In particular, E„ is analytic in Uq. 

Proof. Since E„ is analytic in a neighborhood of with an at most algebraic singularity 
there, it is sufficient to analyze its limit as z — > from the upper half plane. By (44) and 

hm D{z,E,)f{z)-^= hm c^+i ^o^i^^-^ lo^ifi^)) = cyH-\ 

z—>0 z — '0 

lmz>0 Im2>0 



On the other hand, by (43) and (48) , 



lim D (z, w) W (z)-^ = c-i/2gfi-(0) ^ ^-1/2 (i^j^ vr - ih{0)] 



lmz>0 



with <I> given by (|42| and h defined in 1^. 
Summarizing, 



lim D(z,Wcr^W(z) f(zf = —e 

z-»0 c 
lmz>0 



-j$(0) 



By (45) and (69), if Imz > 0, 



with 



En{z) = D2A{z)mn{zr', 



D{z,Wc) 



(70) 
(71) 
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Gathering the hmits computed above, and using that 

hm a{z) = e^^/^ 



lmz>0 



and the definition of 7/„, the statement follows. □ 



Therefore, by construction the matrix- valued function Pq given by (50) satisfies conditions 
(Pol)-(Po4). Moreover, it is easy to check that 

detPo (z) = 1 for every z E C/o\S. 
2.7 Final transformation 

Recall that matrices N and P^, C £ {"1)0,1} have det = 1 in their domains of definition. 
We may define 

zGC\{SU[/_iUt/oU[/i}; 
zGf/c\S, CG {-1,0,1}. 

R is analytic in C\ {T, U dU-i U dUo U dUi}. In fact, since N matches the jump of S on 
( — 1,1), and matches the jumps of S within C/^, £ { — 1,0,1}, we conclude that R is 
analytic in the complement to the contours Tir depicted in Fig. |4] with additional possible 
singularities at {—1,0, 1}. But taking into account (S5) and the local behavior of Pi;^ at these 
points (see (P^4)), we conclude that these singularities are removable. 





Figure 4: Contours E/j. 
Now we compute the jumps of R. For the sake of brevity, we denote 

Sfl \ {dU.i U dUo U dUi) . 



Then by (S2) and ([72]), for z G SfJ**, 

R+ (z) = R_ (z) N (z) 



-2n 



{z) 



(73) 
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On the other hand, for dUj {j £ {—1,0,1}) oriented clockwise, we have that R_|_ (z) 
S+ (z) (z) and R_ (z) = S- (z) PJ^ (z). Hence, 



R+ (z) = R_ (z) P,- (z) N-^ (z) , zedUj, jG {-1,0,1}. 



(74) 



Summarizing, R defined in ([72| is analytic in C\S^, satisfies the jump relations (|73|)-(74) 
on Sr, and has the following behavior as z ^ oo: 



R(z) =1 + 



By (74) and (Po3), as n ^ oo, 

R+(z) = R_ (z) (l + O 



uniformly on dU-i U dUo U dUi. 



(75) 



On the other hand, there exists a constant < g < 1 such that ly? (z;)| ^ < q < 1 uniformly 
on S^*. Since N does not depend on n, we conclude from (73) that as n ^ oo. 



R+(z) = R_ (z) (I + O (g2")) uniformly on 



(76) 



Motivated by (|73|) -(|76j) we define 

1 



A(.) 



N(s) 



Wcis) V(s) 



-2n 



[Pc(s) N-i(s)-I, 



(s) - I, for s € S^"*; 

forsGa?7o iG {-1,0,1}, 



so that R+(z) = R_(z)(I + A(z)), z G S^j. Following |i l4,. Section 8] we can show that A 
has an asymptotic expansion in powers of 1/n of the form 



^ (•5) ~ — ^"^r — ' as n ^ oo, uniformly for s G S^j. 



(77) 



fe=i 



By (76), for A; G N, 



Afc(s) = 0, forsGS?j"*. 



(78) 



Furthermore, by [14" formulas (8.5)-(8.6)], 
(a, — 1) 



Afc(s) 



N(s) e^^c2W{s) 



^^(a^ + ifc-l) -(A:-i). ^ 



2'=[log(/9(s)]'= 

2 c2W{s) N (s), for iblms > and s G at/i 
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and 



-J- iTT0 _1 

where (a, 0) == 1, 

(a, k) 



-0-3 



for iblms > and s G dU^i, 



{4a^ - l){4a^ - 9) • • • (4a^ - (2fc - 1)^) 



Each Afc on the small contours encircling ±1 is independent of n and possesses a meromorphic 
continuation to U-i and Ui with the only pole at ±1 of order at most [(k + l)/2]. However, 
unlike in the case analyzed in [Ij], the existence of a jump in the weight is revealed through 
the contribution of the local parametrix Pq, and hence, each is in general not independent 
on n, although uniformly bounded in n. 

So, it remains to determine A.^ on OUq. Here we calculate explicitly only the first term, 

Ai. 



Using (|45j), (|48j), (|50|), (|63j), (|67|) and (|69j), we obtain 
A 



A (s) = E„(s) 
Let us define 



-A 



-TA 



nf{s) V-l/n A 



+ 



1 



E„^(s), s e dUo, n ^ 00. 



Ai(.)^=^-^E4s) 



-A -Tx 
-1/tx a 



E-\s), s G dUo. 



(79) 



Using that by ( 69 ) , 



E„(s) =F(s) (i"n^ 



0"3 



F(s) f"c 



0-3 



where 



F(.) 



if Ims > 0; 



,0-3 



1' 
-1 0, 



fisf"', iflms<0, 



we conclude that for s £ BUq, 
Ai(s,n) ^ 



F(s) 



f{s) '''' V(-l)"+^n-2Vr, 



A 



(80) 



which is uniformly bounded in n, so that Ai in (79)-(80) is genuinely the first coefficient in 
the expansion (77). 
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Similar analysis can be performed for A./, (^n), k > 2, taking higher order terms in the 
expansion of ^ in ( 62 ) . 



The explicit expression (80 ) and the local behavior of / show that Ai (s, n) has an analytic 
continuation to Uq except for the origin, where it has a simple pole. Again, a similar conclusion 
is valid for other A^ (s, n), except that now the pole is of order k. 



As in [5j Theorem 7.10] we obtain from (77) that 



R(j) {z, n) 



as n — > oo, 



(81) 



uniformly for z G C\ {dU-i U dU^ U dUi} where each R^-') (z) is analytic, uniformly bounded 
in n, and 



O I - I as z — > oo. 



This is a bona fide asymptotic expansion near infinity, since 



V/ G N 3C > : Izl > 2 



R(z) 



R(J) (z,n) 



< 



C 



\z\ n 



for any matrix norm ||-||. The proof is based on the integral representation for R, 



R(z) =1 + 



1 
27ri 



R_(s)A(s,n) 



ds, 2; G C \ Si 



(see [5]); although in our case the coefficients Afc and R^ in (77) and (81) depend on n, 
their uniform boundedness allows to follow the steps of the proof of Lemma 8.3 in [14 . In 
particular, expanding the jump relation R-|- = R_ (I + A) up to order 1/n we find that 

R+^ (s, n) - rL^^ (s, n) = Ai (s, n) , for s G dU-i U dUo U dUi. 

Since R^^) is analytic in the complement of dU^i U dUo U dUi (see (78)) and vanishes at 
infinity, by the Sokhotskii-Plemelj formulas. 



rW (z,n) 



1 
2Td 



Ai (s,n) 



ds. 



I dU-iUdUoUdUi S - ^ 

Recall that Ai can be extended analytically inside C/j's with simple poles at ±1 and 0; let us 
denote by A^^'^ (n), B^^^ (n) and C^^^ (n) the residue of Ai(-,n) at 1, —1 and 0, respectively. 
Then residue calculus gives 



rW (z,n) 



f (n) ^ 5(1) (n) ^ C7(i) (n) 
z — 1 z + 1 z 



for z G C\ {[/_! U [/o U C/i} ; 



AW{n) 5(1) (n) C7«(n) ^ , , , ,,,,,, 

V + + ^ - ^1 -z, " , for z G U-i U f/o U C/i. 

z — 1 z + 1 z 



(82) 
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Residues A^^^ (n) and B^^^ (n) are in fact independent of n; they have been determined in 
[H Section 8]: 



(83) 



(notice however an extra factor ^/c in the constant D^o with respect to jl4j). The value of 
the remaining residue C^-^^ (n) is given in the following 



Proposition 17 Coefficient C^^'^ (n) in (82) is given by 



rW (r,\ n^'^ ( -cos^" A-isin6„\ 3 

27r °°VA + isin0„ cos^„ ' ° 



where On is defined in ([s]). 



Proof. Taking into account (64) and (79) we conclude that 



c<»w = |e„(o)(_-;^^ 7)e„-.(0). 



By Proposition 16 



oo 



With the notation ([7| and choosing ^ € M such that e*'' = t;^, we get 

1 l^ „ir,„<73 -A -e*^\ ^-i^„<x3 A -l^ _ o - cos(27?„ + q) A - i sin(2r/„ + q] 



-1 i; V-e~'^ ^ J Vl 1 y V^ + ^sin(2?7„ + ?) cos(2r?„ + ?) 

It remains to observe that 2r]n + = 9n, and this settles the proof. □ 

3 Asymptotic analysis. Proof of Theorems 

Unraveling the transformations Y^T^S^Rwe can obtain an expression for Y. We 
specify the following domains (see Fig. [5]): 

• De is the unbounded component of C\T,fi; 
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P- correspond to the portion of the inner domain exterior to C/^, C G {— 1)0,1}, lying 
in the upper (resp., lower) half-plane; 



is the subset of in the outer domain and upper (resp., lower) half plane; 
V^- is the subset of f/^ in the inner domain and upper (resp., lower) half plane. 



From (311, (35|, and (72), 



Y{z,n) 



with c e {-1,0,1}. 



1 0^ 



ip[z 



(84) 




Figure 5: Domains for Y 



Next, using the asymptotic expression for R derived above, we obtain information about 
the behavior of Y in different domains of the plane. 

3.1 Asymptotics for the monic orthogonal polynomials on C\[— 1, 1]. Proof 
of Theorem [H 



If i^T is a compact subset of V^, then by (81 ) and (84), 



Y{z, n) = 2-"'^3i^(z)N(z)^"'^3 (z) , z£ K. 



(85) 



Since P„ {z) = Yn (z, n), we get by (|45|-(|47]) that 

2"Pn {z) ^ 
(^(z)" D{z,Wc) 
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Au{z)m{z), 



with 



(86) 



By (47) and (81), uniformly on K, 

2"P„ {z) _ ^ [zf/^ 



{zr D{z,w,)^^^2_ifl^ 



1 + 



Tin [z) 



n 



+ 4 



as n — > oo, 



with 



(87) 



Taking into account the expression for R(i) in ( |82p , as well as ( |83[ ) and Proposition 
get that 



17 



we 



(1) l-4a^ 4/?2 - 1 _ log(c) cos(e„) 



R['^{z) = iD 



16(z-l) 16(z + l) 27rz 

4a2 -1 4/32-1 log(c) log(c) /vr - sin(e„) 



^16(z-l) 16(2 + 1) 2tt 
The trivial identity ^'^{z) + 1 = 2zip(z) yields 



1 ± 



z±l\ ifiz)] ip{z)±V 
and we conclude that in ii' C T^e, TZ^i {z) = T~in {z)^ with Tin defined in ([o]). 



(89) 



3.2 Asymptotics of the recurrence coefficients 

Recall that monic polynomials P„ satisfy the recurrence relation 

Pn+l{x) = {x- hn) Pn{x) - a^P„_i(x), 71 = 0, 1, ... , 

with P_i {x) = and a„ > 0. From [10^ (see also [8J) it follows that the coefficients can be 



found directly in terms of the matrix Y in ( 85 ) : 
a\= lim z'^Yi2 {z, n) Y 21 {z,n 



lim 



Dl 



+ zRi2 {z, n)^ ^zR2i {z, n) + 



(90) 



lim (z - Yii (z, n + 1) Y22 (2, n)) = lim z (1 - Rn {z, n + l) R22 (z, n)) . (91) 

+00 2^00 
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We may take limits in the asymptotic expansion (81); additionally to ( |88| ) we have that 



4a2 _ 1 



Vl6(z - 1) 16(z + 1) ^vr 
4a^ - 1 4/3^ - 1 log(c) cos 
lQ{z - 1) ~ 16(z + i; 



4/3^-1 ^ log(c) log(c)/7r + sin(g„: 



+ 



2ttz 



Thus, 



1 log c 



2im 



sin(^„) + O ( ^ ) , n ^ cx), 



which proves (10). Analogously, 



K = 



_ logC COs(6'n+l) - COs(6'n) q ( ^ 



27r 



n 



n — > oo. 



By ^ and (§, 



9„+i-0„ = 7r + 2i^ log( 1 + 
vr 



n 



so that 



logc cos(g„ + 2^ iog(n±i)) + cos(g„ 
' 27r 
log c 



n 



+ 



27rn 



2 cos (6l„) + O 



+ 



1 



2 ' 



which proves ( 11 ). 



In [T7] A. Magnus conjectured that for the weight 



w (x) = (1 — x)" (1 + x)^ |xo — xl''' 



(92) 



I i?, for X S [— 1, xo) , 
[A, for X e [xo, 1] , 

with A and B > Q and a, /3 and 7 > — 1, and xq G (—1,1), the recurrence coefficients of the 
corresponding orthogonal polynomials exhibit the following behavior n ^ 00: 

1 M 



2 n 
2M 



cos (2nto — 2// log (4n sin to) — + o(l/n) , 



bn = cos y{2n + l)to - 2/ilog (4nsinto) - +o(l/n) 



(93) 
(94) 



where 



IB 1/7^ 

Xo = COs(to), < to < TT, ^=^^°g^' ^ = 2 Y X 



7 



vr - (a + /3 + 7) to - 2 arg r g + i/i 



7 

arg ( - + z/i 



35 



log c 



Taking B = 1, A = c^, j = 0, and xq = (to = 7r/2), we get 
I logc|/(27r), and 

~ = i^-^P)^ - 2 arg r (-A) - arg (-A) 

= ^"~/^'' +2argr(A) + | sgn(log(c)). 

Replacing these expressions in ( 93 ) and using the definition in ([7| we obtain 

1 I log c| 



iA, M 



2 27rn 
1 log c 



cos - ^ sgn(log(c))^ +o(l/n) 



and 



2 27rn 
I log c 



sin(e„) +o(l/n) 



vr 



sin 6'n - - sgn(log(c)) + o(l/n) 
irn \ 2 



log c 



vrn 



cos (6'n) +o(l/n) . 



(95) 



(96) 



Comparing these expressions with (|10[)-( 11 ) we see that Magnus' conjecture is valid for 7 = 0; 
moreover, we have shown that in this situation we can replace the error term o(l/n) by a 
more precise 0{l/n'^). 



3.3 Asymptotics for the the leading coefRcient kn 

1 



By m, 



and with ( 84 1 , 
1 



lim 
27ri z^oo 



^ . lim z-"Y2i(z,n + l) 
27rz z^oo 



n+l 



(zRai (z, n + l) Nil {z) + ZR22 (2, n + l) N21 {z)) 



Taking into account (N3) , ( 47 ) and ( 92 ) , we see that 



lim 2;R2i (z, n + 1) = ^„ 
2^00 nUt^, 



2a2 + 2/?2 - 1 log(c) /log(c 



8 



+ 



27r 



vr 



+ sin(^„+i) +0 



n^ 



lim -zN2i (2;) 

2— >oo 

and 



21)2 



2a2 + 2/32-1 log(c) /log(c) 



+ 



vr 



vr 



+ sin(6'„+i 



1 / 1 
- + 



n 



n^ 



and this proves Theorem [5} 
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3.4 Asymptotics for the monic orthogonal polynomials in Uq and on {—6, 6) . 



By analyticity of PnS, it is sufficient to consider z G Vq- and Rez > 0. Using formulas (50) 
(l66|) and (IMl) we get 



Y {z, n) = l-'^'m (z) En (z) * (n/ (z)) W (z)"'"^ f (z) 



^09-2" 1 



^{zrK (97) 



We are interested in the ffi'st column of Y, which is obtained multiplying the r.h.s. of (97) 
from the right by the column vector (1,0)"^. Observe that 



Wizy^^iz) 



1 

1 ,„^»^-2n 



ifiz)- 



ip{z) 



1/W{z) 
W{z)/w,{z) 



1 



W{z) \l/c 



1 



where we have taken into account the definition of W in "Dq-. Thus, 



W{z)Y {z, n) Q = (^) E„ (^) (^)) 



(98) 



Notice that V^^ is mapped by / onto the sector denoted by ® in Figure jsj and vector (1, 1/c)-^ 
corresponds to the ffi'st column of the jump matrix J2 in (54). Taking into account (^'2) we 



the ffist column of ^ in (57): 



conclude that the product of the last two matrices in the right hand side of ( 98 ) is equal to 

(99) 



*(n/(z)) 



1 

1/c 



r(l-A)G(A;n/(z)) 
r(l + A)G(l + A;n/(z)); ' 



By (70), 



W{z)Y{z,n) 



2-"'^«R(z)DS?A(z)m„(z) 



r(l-A)G(A;n/(z)) 
r(l + A)G(l + A;n/(z)) 



(100) 



with A and m„ defined in (46) and (71), respectively. Taking into account formulas (47), we 
conclude that 



2"VF(z)P„(z) = 



Lp{z) 



1/2 



2 (z2 _ 1 



X {9l(z) mn{z)T (1 - A) G (A; nf (z)) + 9^(z) mn{zy^T (1 + A) G (1 + A; nf (z))} , 



where we have used notation (86) and 



5H(z)^='Rn(^);^ + ^Ri2(^). 
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Inserting again (81) we obtain the asymptotic expansion valid uniformly on compact 



11 Di, V / 12 



subsets of C/q. Using the function IZn defined in (87) and introducing 

we rewrite this identity for P„ as 

2^Pn{z)W{z) =D^Au{z) 



n 



1 + - nn{z) + o{\\\ mn{z)T (1 - A) G (A; nf {z)) 



l + iK„W+ofi,^^r(l + A)G(l + A;„/(.)) 



(101) 



(p{z) \ n \n? J J mn{z) 

Let us simplify this expression for the case when z is on the real line. Taking the limit 



z ^ X £ (—5, 5) from the upper half plane, we get by (u|, (p5|, (68) and Lemma 14 



= e*(^(^)+'?"\ for (-5,5), 



with p{x) given in (13), so that on [—5,5), 

mn{x) = mn{x)~^. 

Additionally, we have A = —A and for x G M, on account of formulas (6.1.23) and (13.1.27) 
from [1], respectively, 

r(l + A) =r(l-A), and G (1 + A; ix) = iFi (-A; 1; -ix) e^^/^ ^ G^A^. (102) 
Finally, on [—5,5), 

V9+(x)l/2 / i^_^^x)-^/2 \ ^-»arcsin(x)/2 



(^ll(^))- 



2(x2 - 1)^^ 



^{x^-i)T, 



{Ai2{x))_ 



^(1 _ a;2)l/4 



and by (80), (^82 , (83) and Proposition 17 



(R(i)),, (x) = (r(i) (x))^^ , (R(i))^2 (^) = - {^^'^ 



X 



12 



Gathering all this information in (101) we conclude that locally uniformly on {—5,5), as 
n — > oo. 



2'^P„(x)Ty(x) 



72 Do 



X Re 



(1 -x2)l/4 

TZnix 



1 + 



+ I ^ I I e 
n \ 



arcsin(x)/2^^(^)p (1 - A) G (A; uf (x)) 
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with TZn defined in (87). Observe however that now the exphcit expression for TZn differs 
from Tin defined in ^hy ([82]), in 



nn{z)=nn- (Ai(z,n)) 



'11 



(Ai (z,n))i2 



with Ai {z,n) given in (79). Using (47), (70), and (71), we get 



(Ai(z,n))ii 
(Ai(z,n)) 



f{z)^izf-l 

\2 



A 1 + 



/(^)v'(z)'-i 



y?(2)V viz) 

Txmn{z) - 



T\mn [zf H ^— — 2 

rAm„ (z) 



1 



99 (z 



TxrUniz) Lp{zy 



Hence, by (89) we have that 



nn {z) = Hn (z) + 



f{z) 



A + 



(p{z) 



For further simplification of our formula we may take into account that by formula 
(6.1.29)], for 1, 



.log c 



r (1 - A) = r (1 + A) = Ar (A) = -i^^ |r (a)| e 

Y^log c sinh (log c) 



iargr(A) 

iT(c)e-*^'^s^(^), 



with T(c) given in (14), and we obtain that for x € (—5, 5), 

r(l - A) mn{x) = -iT(c)e-^^''s^(^)e*(^(^)+'?") = -iT(c)e* 
Analogously, 



so that for x G (—6, 6), 



A 



Summarizing, 
2"P„(x)l^(x) 



A H —Txm„ 



ip{z) 



, a;G(-5,(5), 
z log c / log c 



27rarcsin(x) \ vr 



+ e" 



i (2/9(x)+6,i+arccos(x)) 



X Re 



(1 -x2)l/4 

7?.„(x) 



1 + 



+ 



n 
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which proves Theorem [6| 

Furthermore, with the appropriate rescahng and taking into account the local behavior 
of the terms in the right hand side of the asymptotic expression for Pn we easily get the 
assertion of Corollary |8j 

In order to prove Proposition [o] we rewrite (15) as 



7TX 



n 



Poo T(c) 



in— 1 



lFi (A; l;27rix)\ Im 



(103) 



where (5 is the function introduced in (Hi) of Proposition 12 corresponding to a = log(c)/7r. 
Let us consider here only the case c > 1 (the other case can be easily reduced to c > 1 by a 
change of variables x i-^ —x). Then & is strictly increasing in M. If we denote by 



the solutions of 



<cit<---<ci"^<o<ci")<...<Cl"^<.., 



— 0(27rx)^-^ mod(Z), 



(104) 



then by (103), 



lim 



(105) 



where we have used notation (16). Since S5(0) = 0, we have that (^q""^ is given by 

1 



2tt 



<3{2ttx) 



27r 



where {•} is the fractional part of the number, which by strict monotonicity of (5 shows that 

+ k, k£Z. (106) 



m 



In particular. 



[A:,A. + l)9i-©(27rC^") 



27r 



27rCi"^ -2arg(iFi (A;l;27riC, 



.(n) 



where we have used (ii) of Proposition 12 Hence, 



< Ct^ < 1 and Cl."\ < <k + l, k£Z. 

By compactness and diagonal argument, we can always select a subsequence A C N such that 
the following limits exist: 

limCl"^ = a, keZ. 
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By (p6|), 



1 



& 27rC 



(n) 



27r 



and taking limits we conclude that 

Ck - Ck-1 = 1 + ^ ( arg (iFi (A; 1; liriCk)) - arg (iFi (A; 1; 27riCfc-i)) ) • 



(107) 



Let k € N; since arg(iFi (A; l;27riCfc)) is strictly decreasing in [0,+oo), the second term in 
the right hand side of (107) is < 0, so that we conclude that 

< a - a-i < 1, ken. 



By (105), we obtain that 

< liminf - (xj^"^ - x^^\) < limsup - (x^"^ - x^^}^ < 1, ken. 



In the same vein, since arg (iFi (A; 1; 2TTi(k)) is strictly increasing in (— oo, 0), by (107), 

Cfc - Cfc-i > 1, -ke n, 

so that 



lim inf — ( x 

n TT 



4-^ >1, 



-k G N. 



(n) 

Furthermore, observe that for c 7^ 1, the accumulation points of the sequence Q is dense 
in the interval ©-^([0, 27r]). Indeed, by ((tJ) and 



27r 



n ^ log c 



logn + u, V 



log c 



log4 + 



(3 -a h{0) 1 



+ 



vr vr 



argr(A). 



Since c 7^ 1, we can always take b e {2,3} such that (log c) (log 6) /vr^ ^ Q (indeed, otherwise 
we would have that log 3/ log 2 is rational, which is obviously impossible). For such a b, with 
n = 2 6™, m G N, equation (106) is rewritten as 

1 - ^- f„.l°gc^_^ , , logc 



27r 



6 2<^ 



m- 



log 6 + f + 



log2L 



By Kronecker-Weyl theorem (see, e.g. [H Chapter III]), the sequence 



m- 



log c 



log 2 + u + 



log c 



log 2 



is dense in (0, 1), and it remains to use the strict monotonicity of <3. This finishes the proof 
of Proposition l9| 



41 



Remark 18 Obviously, if c = 1, A = 0, iFi (A; 1; Iniy) = 1, and we obtain the clock behavior 



via (105) and (107). 



Regarding Remark 10 it is easy to check numerically that in general 

Ck — Ck-1 / Cfc+1 — Ck, k £ Z, 

which shows that the quasi-clock behavior fails too. Nevertheless, for function y{x) introduced 



in Proposition 12 we obtain applying (28) that y'{x) = 0(1/1x1), oo, from where the 



"clock behavior in the limit" (18) follows by (107). 



3.5 Christoffel-Darboux Kernel 



Using the Christoffel-Darboux formula \21\ Section 3.2], we can write the kernel (19) as 

kn-l Pn (x) Pn-l (y) - Pn (v) Pn-1 (x) 



Kn{x,y) 



k 



kn x-y 

2 Pn (X) Pn-1 (y) - Pn iv) Pn-1 (x) 



n-1' 



X-y 



xy^y; 



using ( 30 ) and the fact that det Y = 1 , we obtain 

-1 Yii (x, n) Y21 (y, n) - Yn {y, n) Y21 (x, n) 



Kn {x, y) 



2ni 
1 1 



x-y 



(0, l)Y-i(y,n)Y(x,n) Q . 



27ri x-y 

By analyticity, it is obviously sufficient to compute Kn when x,y £ (0,6). From (100) 
'1 



(108) 



W{x)Y+ (x,n) 



2-""3R(x)L»SA+(x) rUnixY 



r(l-A)G(A;n/(x)) 
r(l + A)G(l + A;n/ (x)) 



(109) 



On the other hand, since det Y = 1, we have by (97), (69), (45), i71v, 



(y, n)-' = if {yr^--^ {_i}^-2n ?) ^ {vr' W {yr [* (n/ {y)T^ 
X m„(2/)--3A+(y)-iD^'^^^R(y)-i2"'^^ 



The matrix ^ built in (56)-(|6lj) also satisfies det = 1, so that considerations that lead us 
to ( [Tool ) show that 



(0, 1) W(y)Y(y, 



n 



X m. 



-r (1 + A) G (1 + A; n/ (y)) , T (1 - A) G (A; n/ (y))) 
,(y)-'^3A+(y)-iD-3R(y)-i2-3. 



(110) 
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Observe also that locally uniformly for z G (—5, 5), R(z) = I + 0(l/n), which implies that 

R (y) ""^ R (x) = I + ©(l/n), locally uniformly for x,y £ (—(5, 5). 
Gathering ( 109[ ) and (110) in (108) we conclude that for x 7^ y, 



lim - w [ — )W [^) Kn(—^ — ] = — 

n^oo n Vn/ Vn/ Vn n/ 27rz x — y 

(-r (1 + A) G (1 + A; 2„,) , r (1 - A) G (A; 2„,)) L (Vl^j'^i^ '^^.'^^t; 



By formula (6.1.31) of |T and the definition of A in (39), 

logc 2c log c 



r(i + A)r(i-A) 



sinh (log c) 



c2-l 



Since 



TTXX 

n J 



lim ( — ) = ^/ch{0), 



(23). 



we obtain (21 )-(22). Taking into account (102) we can easily rewrite this formula in the form 



Finally, the confluent form of the kernel in (22 ) is obtained from the expression for x ^ y 
by taking limit y ^ x. 



4 Properties of the confluent hypergeometric function 



We prove finally the properties of iFi and related functions summarized in Proposition 12 



(i) Let Kn be the Christoffel-Darboux kernel defined in (19). Then for 2; G C, 

n-l 



Kn{z,z) =Y,\Pk {z)\^>0. 



k=0 



This property is obviously inherited in the limit (21), which implies that 



Koo{z,z)>0, zGC. 
On the other hand, from formula (13.1.27) of [1 it follows that 



i-Fi {ia; 1; iz) = iFi {—ia; 1; —iz) = e {ia + 1; 1; iz) 



G {ia; iz) = G {ia + 1; iz) 



[III) 
[112) 
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Thus, iov z e C\R, X = ia, a £ R \ {0}, 

1 logc G{X;2TTiz)G{l + X;27riz) -G{1 + X;2TTiz)G{X;2TTiz) 



Koo {z, z) 



27r/i(0)c2-l Imz 
1 logc \G{X;2mz) ^ - |G'(A;27riz) p 
27r/i(0) c2 - 1 h^z 



> 0, 



which yields that 



\G{X]iz)\>\G{X]iz)\, lmz>0. 



(113) 
(114) 



Assume that for C G C+, G(A;iC) = 0; then by (112) and (114), 



Hence, 



G(A;0 = G(m + l;iC) = 0. 



iFi (ia; 1; iC) = i^^i (1 + ia; 1; iQ = 0. 



By induction and recurrence relation (13.4.1) in |1J we conclude that every iFi {ia + n; 1; i^), 
with n E Z U {0}, vanishes. But this is impossible, as follows from the addition formula 



iFi(A;l;z + C) 



C 

z + C 



(A)„z" 



A oo 
n=0 ^ 



(115) 



(see ISni formula (2.3.4)]). 



Thus, we conclude that fi{z) = G(A, iz) G HB. The assertion for /2 is obtained by means 



of formula ( 111 ). 



(ii) In order to prove (|26|) assume that a G M \ {0} and x G M. Then by (111), 



iFi (ia; 1; ix) = iFi {—ia; 1; —ix) = e ^\Fi {ia + 1; l;ix) 



(116) 



Hence, an assumption that for a; G M, iFi {ia; 1; ix) = implies that iFi {ia + 1; 1; ix) also 
vanishes, and we arrive at a contradiction reasoning as above and using the addition formula 
( fTl5| ). 

Furthermore, the location of the zeros in the corresponding half planes and the inequality 



(27) is a direct consequence of (i). In particular, function 

iFi (1 + ia; 1; iz) 



h{z) 



iFi {ia; l;iz) 



is holomorphic in C^, continuous in C+ = U M, and satisfies \h{z)\ < 1 for z G C+ and 
|/i(2;)| = 1 for z G M. In consequence, by the maximum principle, \h{z)\ < 1 for z G C"*", 
which proves that the inequality in (27) for z in the upper half plane is strict. 
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(in) Due to (ii), y(x) is correctly defined and real-analytic on M, in particular, y' can 
vanish only at a discrete set of points that can accumulate only at infinity. Again by (13.1.27) 

of m, 



iFi (1 +m;l;ix) ^ iFi(-m;l;-ix) ^ ^i{x-2y{x)) 
iFi (za; 1; ix) iFi{ia;l;ix) 
With the straightforward identity (see (13.4.4) in [1 ) 

iFi (1 + ia; 1; ix) — iFi {ia; 1; ix) = (ix) iFi (1 + ia; 2; ix) 



(117) 



we rewrite (117) as 



J(x-2y{x)) 



1 + ix 



iFi (1 + ia; 2; ix) 



iFi {ia; 1; ix) 

Since the real part of the left hand side is < 1 , this implies that 



Im 



L^i (1 + ia; 2; ix) ' 



1^1 [ia; 1; ix) 
On the other hand, by (13.4.8) of [I], 

y'{x) = Im (-^ log (iFi (ia; I; ix) 



> 0, for X > 0, 
< 0, for X < 0. 



-aim 



iFi (1 + ia; 2; ix) 
iFi (ia; 1; ix) 



(118) 



:ii9) 



(120) 



and by inequality ( 119 ), the first part of the statement of (ii) follows. The differential equation 
in (28) is obtained by taking the real part in (118) and using (120). 

In order to prove (iv) we observe that C5 satisfies the following initial value problem: 

x6'(x) = x + 2a(l-cos0(x)), 0(0) = 0. (121) 

For a = the statement is trivial. Assume first that a > 0; then by (ii), we only need to 
prove that 0'(x) > for x < 0. 

Since & is also real analytic, expanding it at x = we readily conclude from (121) that 
0'(O) = 1. Hence, is locally increasing at the origin. Differentiating (121) we obtain that 



X 0"(x) = 1 + 0'(x) (-1 + 2a sin 0(x)) . 



If for X = C < 0, 0'(C) = 0, then by ( p2| ), 

C0"(C) 



1. 



(122) 



(123) 



In particular, 0"(C) < 0, which shows that every critical point of in the negative semi-axis 
is a strict local maximum, which is incompatible with the behavior at the origin. Thus, 0' is 
sign-invariant on (— oo,0), and in consequence, 0'(x) > there. 

Assume now a < 0; again by (ii), we only need to prove that 0'(x) > for x > 0. 
Reasoning as above, if for ( > we have 0'(C) = 0, then we get (123), which shows that 
every critical point of in the positive semi-axis is a strict local minimum, which is again 
incompatible with the behavior at the origin. 
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Remark 19 The proof (i) presented here is an evolution of an idea of Doron Lubinsky 
to look at the ChristofFel-Darboux kernels Kn for each finite n G N and then taking limits 
as n — > cxd; see also [15^. The argument that allows to conclude that iFi does not vanish, 



based on the addition formula (115), and the identity (117) were suggested to us by Adri 



Olde Daalhuis. We gratefully acknowledge these two contributions of our colleagues. 
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